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Abstract 

The paper focuses on a class of light-tailed multivariate probability distributions. These are 
obtained via a transformation of the marginals from a heavy-tailed original distribution. This class 
was introduced in Balkema et al. (2009). As shown there, for the light-tailed meta distribution the 
sample clouds, properly scaled, converge onto a deterministic set. The shape of the limit set gives 
a good description of the relation between extreme observations in different directions. This paper 
investigates how sensitive the limit shape is to changes in the underlying heavy-tailed distribution. 
Copulas fit in well with multivariate extremes. By Galambos's Theorem existence of directional 
derivatives in the upper endpoint of the copula is necessary and sufficient for convergence of the 
multivariate extremes provided the marginal maxima converge. The copula of the max-stable limit 
distribution does not depend on the marginals. So marginals seem to play a subsidiary role in 
multivariate extremes. The theory and examples presented in this paper cast a different light on the 
significance of marginals. For light-tailed meta distributions the asymptotic behaviour is very sensitive 
to perturbations of the underlying heavy-tailed original distribution, it may change drastically even 
when the asymptotic behaviour of the heavy-tailed density is not affected. 
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1 Introduction 

In recent years meta distributions have been used in several applications of multivariate probability the- 
ory, especially in finance. The construction of meta distributions can be illustrated by a simple example. 
Start with a multivariate spherical t distribution and transform its marginals to be Gaussian. The new 
distribution has normal marginals. It is called a meta distribution with normal marginals based on the 
original t distribution. Since the copula of a multivariate distribution is invariant under strictly increas- 
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ing coordinatewise transformations, the original distribution and the meta distribution share their copula 
and hence have the same dependence structure. The light-tailed meta distribution inherits not only the 
dependence properties of the original t distribution, but also the asymptotic dependency. These asymp- 
totic properties are of importance to risk theory. They include rank-based measures of tail dependence, 
and the tail dependence coefficients (see e.g. Chapter 5 in [8]). Vectors with Gaussian densities have 
asymptotically independent components, whatever the correlation of the Gaussian density. The meta 
density with standard Gaussian marginals based on a spherically symmetric Student t density has the 
copula of the t distribution, and the max-stable limit distributions for the coordinatewise maxima also 
have the same copula. The max-stable limit vectors have dependent components. 

However, this recipe for constructing distributions with Gaussian, or more generally, light-tailed 
marginals, based on a heavy-tailed density with a pronounced dependency structure in the limit, has 
to be treated with caution. The limit shape of the sample clouds from the meta distribution is affected 
by perturbations of the original heavy-tailed density, perturbations which are so small that they do not 
affect the multivariate extreme value behaviour. In going from densities with heavy-tailed marginals to 
the meta densities with light-tailed marginals the dependence structure of the max-stable limit distribu- 
tion is preserved by a well-known invariance result in multivariate extreme value theory. In this paper 
we shall give exact conditions on the severity of changes in the original heavy-tailed distribution which 
are allowed if one wants to retain the asymptotic behaviour of the coordinatewise extremes. It will be 
shown that perturbations which are negligible compared to these changes may affect the limit shape of 
the sample clouds of the light-tailed meta distribution. 

Multivariate distribution functions (dfs) have the property that there is a very simple relation between 
the df of the original vector and the df of the coordinatewise maximum of any number of independent 
observations from this distribution. One just raises the df to the given power. This makes dfs and in 
particular copulas ideal tools to handle coordinatewise maxima, and to study their limit behaviour. This 
rather analytic approach sometimes obscures the probabilistic content of the results. The approach via 
densities and probability measures on R"* which is taken in this paper may at first seem clumsy, but it 
has the advantage that there is a close relation to what one observes in the sample clouds. 

Our interest is in extremes. The asymptotic behaviour of sample clouds gives a very intuitive view 
of multivariate extremes. The limit shape of sample clouds, if it exists, describes the relation between 
extreme observations in different directions; it indicates in which directions more severe extremes arc 
likely to occur, and how much more extreme these will be. It has been shown in [2] that sample clouds 
from meta distributions in the standard set-up, see below, can be scaled to converge onto a limit set. 
The boundary of this limit set has a simple explicit analytic description. The limit shape of sample 
clouds from the meta distribution contains no information about the shape of the sample clouds from the 
original distribution. The results of the present paper support this point. 
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The aim of the paper is to investigate stabiUty of the shape of the hmit set under changes in the 
original distribution. We look at changes which do not affect the marginals, or at least their asymptotic 
behaviour. Keeping marginals (asymptotically) unchanged allows us to isolate the role played by the 
copula. We shall examine how much the original and meta distributions in the standard set-up may be 
altered without affecting the asymptotic behaviour of the scaled sample clouds. This shows how robust 
the limit shape is. Then we move on to explore sensitivity. The limit shape of the scaled sample clouds 
from the light-tailed meta distribution is very sensitive to certain slight perturbations of the original 
distribution, perturbations which affect the density in particular regions. For heavy-tailed distributions 
the region around the coordinate planes seems to be most sensitive. 

The present paper is a follow-up to [2] . The latter paper contains a detailed analysis of meta densities 
and gives the motivation and implications of the assumptions in the standard set-up. It presents the 
derivation and analysis of the limit shape of the sample clouds from the light-tailed meta distribution, 
and may be consulted for more details on these subjects. In the present paper Section 2 introduces the 
notation and recalls the relevant definitions and results from [2]. Section 3 is the heart of the paper; 
here we present details of the constructions which demonstrate robustness and sensitivity of the limit 
shape and the asymptotics of sample clouds from meta distributions. Concluding remarks are given in 
Section 4. The appendix contains a few supplementary results, and a summary of notation which the 
reader may find useful when reading the paper. 
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2.1 Definitions and standard set-up 

The following definition describes the formal procedure for constructing meta distributions. 

Definition 1. Let Gi, . . . ,Gd be continuous dfs on K which are strictly increasing on the intervals 
li = {0 < Gi < 1}. Consider a random vector Z in M'' with df F and continuous marginals Fi, 
i = 1, . . . ,d. Define the transformation 

K{xu...,Xd) = {Ki{xi),...,Ka{xa)), Ki{s) = Fr\G,{s)) i = l,...,d. (2.1) 

The df G = F o K is the meta distribution (with marginals Gi) based on the original df F. The random 
vector X is said to be a meta vector for Z (vdth marginals Gi) if 

Z = K{X.). (2.2) 

The coordinatewise map K = Ki ^ ■ ■ ■ ^ K4 which maps x = {xi,. . . , Xd) G I = Ii x ■■ ■ x 1^ into the 
vector z = {Ki{xi), . . . , Kd{xci)) is called the meta transformation. 
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The class of distributions above is too general for our purpose. Hence, we choose to restrict our 
attention to a subclass by imposing more structure on the original distribution and on the marginals of 
the meta distribution. The standard set-up of this paper is the same as in [2]. Recall the basic example 
we started with in the introduction. 

The multivariate t density has a simple structure. It is fully characterized by the shape of its level sets, 
scaled copies of the defining ellipsoid, and by the decay c/r^^'^ of its tails along rays. The constant A > 
denotes the degrees of freedom, d the dimension of the underlying space, and c is a positive constant 
depending on the direction of the ray. In the more general setting of the paper, the tails of the density 
are allowed to decrease as cL{r)/r^^''' for some slowly varying function L and the condition of elliptical 
level sets is replaced by the requirement that the level sets are equal to scaled copies of a fixed bounded 
convex or star-shaped set (a set D is star-shaped if z e .D implies e .D for < t < 1). Due to the 
power decay of the tails, the density / is said to be heavy-tailed. Densities with the above properties 
constitute the class T\. 

Definition 2. The set T\ for A > consists of all positive continuous densities f on M."^ which are 
asym,ptotic to a function of the form /,(n£)(z)) where f-t{r) = L(r)/r'^'^+'') is a continuous decreasing 
function on [0,cx)), L varies slowly, and is the gauge function of the set D. The set D is bounded, 
open and star-shaped. It contains the origin and has a continuous boundary. 

The reader may keep in mind the case where D is a convex symmetric set. In that case the gauge 
function is a norm, and D the unit ball. The normal marginals of the meta density are generalized to 
include densities whose tails are asymptotic to a von Mises function: go{s) e"'^'*^ for s — > oo with 
scale function a = l/tp', where tp is a. C"^ function with a positive derivative such that 

V'(s) 00, a(s)' -^0 s 00. (2.3) 

This condition on the meta marginals ensures that they lie in the maximum domain of attraction of the 
Gumbel limit law exp(— e~^), a; > 0; see e.g. Proposition 1.4 in [10]. In this case we say that the meta 
distribution is light-tailed. 

Definition 3. In the standard set-up, the density f lies in T\ for some A > 0, and go is continuous, 
positive, symmetric, and asymptotic to a von Mises function e^^ . We assume that tp varies regularly in 
infinity with exponent > 0. The density g is the meta density with marginals go, based on f. 

2.2 Convergence of sample clouds 

An n-point sample cloud is the point process consisting of the first n points of a sequence of independent 
observations from a given distribution, after proper scaling. We write 

Ar„ = {Zi/a„,...,Z„/a„} (2.4) 
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where Zi, Z2, . . . are independent observations from the given probability distribution on R'^, and a„ are 
positive scahng constants. It is custom to write Nn{A) for the number of the points of the sample cloud 
that fall into the set A. 

In this section, we discuss the asymptotic behaviour of sample clouds from the original density / and 
from the associated meta density g in the standard set-up. The difference in the asymptotic behaviour is 
striking: sample clouds from a heavy-tailed density / converge in distribution to a Poisson point process 
on M*^ \ {0} whereas sample clouds from a light-tailed meta density g tend to have a clearly defined 
boundary. They converge onto a deterministic set. 



2.2.1 Convergence for densities in and measures in Vx 

For densities in J^x, A > 0, there is a simple limit relation: 

fir-n-Wn) _^ w„ ^ w / 0, r„ ^ 00, (2.5) 

J*{rn) 

where 

/i(w) = l/nr,(w)^+'^ = r]{u)/r^+'^, r = ||w||2 > 0, to = w/r. (2.6) 

Convergence is uniform and on the complement of centered balls. If Zi,Z2,... are independent 
observations from the density / then the sample clouds Nn in (2.4) converge in distribution to the 
Poisson point process with intensity h weakly on the complement of centered balls for a suitable choice 
of scaling constants o„. 

A probability measure tt on varies regularly with scaling fimction a{t) ^ 00 if there is an infinite 
Radon measure p on R'^ \ {0} such that the finite measures tn scaled by a{t) converge to p vaguely on 
M'' \ {0}. The measure p has the scaling property 

p{rA) = p{A)/r^ r>0 (2.7) 

for all Borel sets A in M'' \ {0}. The constant A > is the exponent of regular variation. If it is positive 
then p gives finite mass to the complement of the open unit ball B and weak convergence holds on the 
complement of centered balls. We shall denote the set of all probability measures which vary regularly 
with exponent A > by Vx- In particular tt G Vx if tt has a continuous density in J-x- As above for 
independent observations Zi,Z2,... from a distribution tt G Vx the scaled sample clouds Nn in (2.4) 
(with scaling constants a(n)) converge in distribution to a Poisson point process with mean measure p 
weakly on the complement of centered balls (since the mean measures converge, see e.g. Proposition 3.21 
in [10] or Theorem 11.2.V in [5]). 
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2.2.2 Convergence for meta densities in the standard set-up 

Sample clouds from light-tailed meta densities in the standard set-up, under suitable scaling, converge 
onto a deterministic set, referred to as the limit set, in the sense of the following definition. 

Definition 4. For a compact set E in M**, the finite point processes Nn converge onto E if for open sets 
U containing E, the probability of a point outside U vanishes, F{Nn{U^) > 0} — > 0, and if 

P{Ar„(p + eB)>m}^l m>l,e>0,pGE. 

We now recall Theorem 2.6 of [2], which characterizes the shape of the limit set for meta distributions 
in the standard set-up. 

Theorem 2.1. Let f, g and go e~'^ satisfy the assumptions of the standard set-up. Define 

E := Ex,e = {u e M'^ | + • • • + \ud\' + A > (A + d)|lu||^}. (2.8) 

If ^(.^n) ~ logn, then for the sequence of independent observations Xi,X2, . . . from the meta density g, 
the sample clouds Mn = {Xi/r„, . . . , X„/r„} converge onto E. 

2.3 Further notation and conventions 

In order to ease the exposition, wc introduce some additional assumptions and notation. 

All univariate dfs are assumed to be continuous and strictly increasing. The dfs Fq and Fq on R are 
tail asymptotic if 

Fo{-t)/Fo{-t) ^1 (1 - Fo(t))/(l - Fo{t)) ^1 t ^ oo. 

The sample clouds from a heavy-tailed df F converge to the point process JV if iV is a Poisson point 
process on M*^ \ {0}, and if the sample clouds converge to N in distribution weakly on the complement of 
centered balls, where the scaling constants c„ satisfy 1 — -Fb(cn) ~ ^/n. Two heavy-tailed dfs F* and F** 
have the same asymptotics if the marginals are tail asymptotic and if the sample clouds converge to the 
same point process. The light-tailed dfs G* and G** have the same asymptotics if the marginals are tail 
asymptotic and the sample clouds converge onto the same compact set E*, with scale factors 6„ which 
satisfy — log(l — Go{bn)) ~ logn. One may replace a scaling sequence a„ by a sequence asymptotic to 
an without affecting the limit. The scaling of sample clouds is determined up to asymptotic equality by 
the marginals. Tail asymptotic marginals yield asymptotic scalings. 

3 Results 

We now turn our attention to the main issues of the paper. The aim is to investigate how much the dfs F 
and G = F o K in the standard set-up may be altered without affecting the asymptotic behaviour of the 
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scaled sample clouds. For simplicity we assume here that the marginal densities of F also all are equal to 
a positive continuous symmetric density /q. It follows that the components of the meta transformation 
K are equal: 

K ■.x^z={Ko{xi),...,Koixd)) Ko = Fo^oGo Ko{-t) = -Koit). (3.1) 

Let F* be a df with marginal densities fo- Then G* = F* o K is the meta distribution based on the df 
F* with marginals go- One can pose the following questions: 

(Q.l) If the scaled sample clouds from G* and from G converge onto the same set E, do the scaled sample 
clouds from F* converge to the same point process A'' as those from F? 

(Q.2) If the scaled sample clouds from F* and from F converge to the same Poisson point process A'', do 
the scaled sample clouds from G* converge onto the same set E as those from G? 

The answer to the corresponding questions for coordinatewise maxima and their exponent measures 
(if we also allow translations) is "Yes" . Here, for sample clouds and their limit shape, the answer to both 
questions is "No" . This section contains some counterexamples which will be worked out further in the 
next two sections. 

If we replace / by a weakly asymptotic density /* x /, the asymptotic behaviour of sample clouds 
from g* is not affected, since g* >^ g (see Lemma 3.2 below), but the scaled sample clouds from /* 
obviously need not converge. What if they do? 

Example 1. Let Z have a spherical Student t density /(z) = /»(||z||) with marginals /o and limit fimction 
/i(w) l/||w||^+'' with marginals c/t^^^ . The vector with components Z, . . . , aJZ where ai, . . . , 
are independent unit vectors will have the same marginals /o but density /*(n£;(z)) with elliptic level 
sets, which arc spherical only if the vectors a.j arc orthogonal. There are many star-shaped sets D for 
which /i(w) — l/ri£)(w)^+'' has marginals c/t^^^ as above. Probability densities / ^ f*{nD) will lie in 
J-\ and have marginals asymptotic to the Student t marginals /o above. All these densities are weakly 
asymptotic to each other, f ^ f . The only difference between the densities is in their copulas. The 
information of the dependency contained in the set D is preserved in the limit of the sample clouds from 
the density /, but lost in the limit shape E of the sample clouds from the meta density g. Surprisingly 
the information on the shape is lost in the step to the meta density g, but the tail exponent A of the 
marginals is still visible in the limit shape E. 

What we want to do is fix the marginals /o and go (which determine the meta transformation K), and 
then vary the copula and check the limit behaviour of the sample clouds (where we impose the condition 
that both converge). We are looking for dfs F* and G* with the properties: 

(P.l) F* has marginal densities /o; 
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(P.2) G* is the meta distribution based on F* with marginal densities go; 
(P.3) the scaled sample clouds from F* converge to a Poisson point process N*; 
(P.4) the scaled sample clouds from G* converge onto a compact set E*. 

Moreover one would like E* to be the set E\,e in (2.8), or N* to have mean measure p* ~ p with 
intensity h in (2.6). So we either choose F* to have the same asymptotics as F, or G* to have the same 
asymptotics as G. Note that the four conditions above have certain implications. The mean measure p* 
of the Poisson point process N* is an excess measure with exponent A, see (2.7), its marginals are equal 
and symmetric with intensity A/jij'^"'""'^ since the marginal densities /o are equal and symmetric and the 
scaling constants c„ ensure that p{wd > 1} = 1- The limit set E* is a subset of the cube C — [— 1, 1]"^ 
and projects onto the interval [—1, 1] in each coordinate, again by our choice of the scaling constants. 

The two sections below describe procedures for altering distributions without changing the marginals 
too much. A block partition is a special kind of partition into coordinate blocks. If the blocks are 
relatively small then the asymptotics of a distribution do not change if one replaces it by one which 
gives the same mass to each block. Block partitions are mapped into block partitions by K. The mass 
is preserved, but the size and shape of the blocks may change drastically. The block partitions provide 
insight in the relation between the asymptotic behaviour of the measures dF* and dG* . In the second 
procedure we replace dF by a probability measure dF, which agrees outside a bounded set with d(F + F°) 
where F° has lighter marginals than F: 

F^i-t) «Fo{-t) 1-F°{t) «l-Fo{t) t^oo, j = l,...,d. (3.2) 

This condition ensures that F and F have the same asymptotics. The two corresponding light-tailed 
meta dfs G and G on x-spacc may have different asymptotics since the scaling constants 6° and 6„ may 
be asymptotic even though G° has lighter tails than G. If this is the case, and the scaled sample clouds 
from G° converge onto a compact set E°, then those from G converge onto the union EUE°, which may 
be larger than E. These two procedures enable us to construct dfs F* with marginal densities /□ and 
meta dfs G* = F* o K with marginal densities qq which exhibit unexpected behaviour: 

(Ex.1) G* and G have the same asymptotics, but the scaled sample clouds from F* converge to a Poisson 
point process which lives on the diagonal. (See Theorem 3.7 and 3.8, and Example 2.) 

(Ex.2) The scaled sample clouds from G* converge onto E* = Al^ Eqq, where i^oo is the diagonal cross 

i^oo = H I < r < 1, (5 e {-1, l}-^}, (3.3) 

and A C [—1, l]"* any compact star-shaped set with continuous boundary. The dfs F* and F have 
the same asymptotics. The density /* is asymptotic to / on every ray which does not lie in a 
coordinate plane. (See Example 2.) 
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What does the copula say about the asymptotics? Everything, since it determines the df if the 
marginals are given; nothing, since the examples above show that there is no relation between the 
asymptotics of F* and the asymptotics of G* even with the prescribed marginals /o and go. One might 
hope that at least the parameters A and 6, determined by the marginals, might be preserved in the 
asymptotics. The point process N* reflects the parameter A in the marginal intensities A/|t|^+^. However, 
E* = E\*fi* may hold for any A* and 0* in (0, oo) by taking A = E\*^0* in the second example above. 

We now start with the technical details. 

The construction procedures discussed above will change an original df F with marginals Fq into a df 
F whose marginals Fj are tail equivalent to i<o. This is no serious obstacle. 

Proposition 3.1. Let the scaled sample clouds from F converge to a point process N , and let the scaled 
sample clouds from G = F o K converge onto the compact set E. If the marginals Fj are continuous and 
strictly increasing and tail asymptotic to Fq then there exists a df F* with marginals Fq such that F* has 
the same asymptotics as F and G* = F* o K has the same asymptotics as G. 

If moreover F has a density f with marginals asymptotic to fo in ±oo, then F* has a density /*, and 
for any vector w with non-zero coordinates and any sequence w„ w and r„ ^ oo there is a sequence 
^ w such that f*{rnWn) ~ firnyf'J- 

Proof Let F* = F o Kp be the meta df based on F with marginals Fq. The components Kpj = Fj o Fq 
are homeomorphisms and satisfy Kpjit) ~ t for \t\ oo. (Here we use that the marginal tails vary 
regularly with exponent —X^O.) It follows that 

||i^Hz)-z||/||z||-0 H-oo. (3.4) 

This ensures that F and F* have the same asymptotics. (For any e > the distance between the scaled 
sample point Z/c„ and Kp{Z)/cn is bounded by e||Z||/c„ for ||Z|| > ec„ and n > n^.) A similar argument 
shows that G = F o K and G* = F* o K, the meta df based on G with marginals Go, have the same 
asymptotics. Here we use that Gj = Fj o Ko is tail asymptotic to Go = -Fb ° Kq since Fj is tail asymptotic 
to Fq. Under the assumptions on the density the Jacobean of Kp is asymptotic to one in the points 
r„w^ and (3.4) gives the limit relation with r„w^ = ii'^^(r„w„). ^ 

In general, the densities /* and / (in the notation of the above proposition) are only weakly asymptotic, 
as in Proposition 1.8 in [2]. The density /* on z-space is related to / in the same way as the density g* 
is related to g. If /* >: / or /* < Cf or /* /, then these relations also holds for g* and g, and vice 
versa. Similarly for the marginals: gj gfo in oo implies fj /o in oo. These results are formalized in 
the lemma below. 

Lemma 3.2. If F* has density f* with marginals f* and G* = F* o K has density g* with marginals 

g*j, then 

5*(x)/3(x) = r(z)//(z) g*j{s)/go{8) = f*j{t)lfo{t) z = K{x), t = Ko{s). 
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Proof The Jacobean drops out in the quotients. % 
3.1 Block pcirtitions 

We introduce partitions of M'^ into bounded Borel sets i3„. In our case the sets are coordinate blocks. 
Since our dfs have continuous marginals the boundaries of the blocks are null sots, and wc shall not 
bother about boundary points, and treat the blocks as closed sets. To construct such a block partition 
start with an increasing sequence of cubes 

SnC =[-Sn,Sn\'^ < Si < S2 < • • • , S„ ^ OO C=[-l,l]'^. 

Subdivide the ring i?„ = Sn+iC \ s„C between two successive cubes into blocks by a symmetric partition 
of the interval [— s„+i, s„+i] with division points ±s„j, j = 1, . . . , nin, with 

This gives a partition of the cube Sn+iC into (2m„ + 1)** blocks of which (2m„ — 1)** form the cube 
SnC. The remaining blocks form the ring The meta transformation K transforms block partitions in 
x-space into block partitions in z-space. A comparison of the original block partition with its transform 
gives a good indication of the way in which the meta transformation distorts space. 

Definition 5. A partition ofU.'^ into Borel sets An is regular if the following conditions hold: 

1) The sets An are hounded and have positive volume \An\ > 0; 

2) Every compact set is covered by a finite number of sets An; 

3) The sets An are relatively small: There exist points p„ e An with norm ||pn|| = r„ > such that 
for any e > 0, ^„ C p„ + er„B, n > n^. 

The block partition introduced above is regular if and only if .s„-|-i ^ s„ and A„/s„ where A„ 
is the maximum of s„i, s„2 — Snij • • • , s„m^ — Snnir,-!- Regular partitions give a simple answer to the 
question: If one replaces / or g by a discrete distribution, how far apart arc the atoms allowed to be if 
one wants to retain the asymptotic behaviour of the sample clouds from the given density? 

Proposition 3.3. Let Ai,A2, ... be a regular partition. Suppose the sample clouds from the probability 
distribution fi scaled by rn converge onto the compact set E. Let jl be a probability measure such that 
jl{An) = /u(^n) for n> no- Then the sample clouds from fb scaled by rn converge onto E. 

Proof Let p & E, and e > 0. Let /x„ denote the mean measure from the scaled sample cloud from fi and 
fin the same for jl. Then y^„(p + {e/2)B) — > oo. Because the sets An are relatively small there exists ni 
such that any set An which intersects the ball r„p + r„e_B with n > ni has diameter less than ernf^- Let 
Un be the union of the sets An which intersect r„p + (r„e/2)_B. Then [/„ C r„p + er„i? and hence 

/"n(P + (e/2)-B) < lln{Un/rn) = fin{Un/rn) < ftn{P + ^B). 
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Similarly jln{U'') — * for any open set U which contains E. ^ 
Remeirk 1. The result also holds if /i(A„) x provided is positive eventually. 

There is an analogous result for regular partitions in z-space. 

Proposition 3.4. Suppose n G T^x{p) with scaling constants Cn- Let Ai,A2,... he a regular partition 
and let tt be a probability measure on M** such that 7r(A„) = 7r(j4„) for n > no. Then w G T^\{p) with 
scaling constants Cn ■ 

Proof Any closed block A C Mf^ \ {0} whose boundary carries no p-mass is contained in an open block 
U with p(U) < p{A) + e. As in the proof of the previous proposition for n > ni there is a union t/„ of 
atoms An such that A c Un/cn C U. ^ 

For excess measures p with a continuous positive density h there is a converse. For tt e with 
limit p there is a regular partition Ai,A2, . . . such that 7r(A„) ^ /(z)dz with / e J^x, and with the 
same scaling constants. See [4], Theorem 16.27. This result vindicates our use of densities in J^x. Not 
every distribution in the domain of h has a density, or even a density in J^x, but every such distribution 
is close to a density in J^x in terms of a regular partition. 

We thus have the following simple situation: Ai,A2,... is a block partition in x-space and Bi = 
K{Ai), B2 = K{A2), . . . the corresponding block partition in z-space. Let tt be a probability measure in 
z-space and fi a probability measure in x-space, linked by K, i.e. tt = K{jl). Then 7ir(^„) = jl{Bn) for 
all n. So 

7t{An) ~ / /(Z)dz ^ fi{Bn) ^ / g(x)dx. (3.5) 

Theorem 3.5. // both partitions are regular and one of the equivalent asymptotic equalities in (3.5) 
holds, then the sample clouds from tt scaled by Cn converge to the Poisson point process with intensity h 
in (2.6), and the sample clouds from fi scaled by rn converge onto the set E = Ex,e in (2.8). 

Proof Combine Proposition 3.3 and 3.4. ^ 

Unfortunately the meta transformation K is very non-linear. Regularity of one block partition does 
not imply regularity of the other block partition. 

We first consider the case when the block partition (A„) in x-space is regular, but (-B„) is not. The 
block partition {An) in x-space is based on a sequence of cubes s„C = [— Sn, Sn]''- Successive cubes are 
of the same size asymptotically, s„+i s„. The cubes tnC in z-space with tn = -ft'o(sn) may grow very 
fast. It is possible that tn « tn+i, as in Proposition 3.6 below. The corresponding partition with blocks 
Bn = K{An) in z-space then certainly is not regular. 

Proposition 3.6. Let e G (0,1). There is a sequence < Si < S2 < ■ ■ ■ such that s„ ^ oo and 
Sn+i ~ Sn, and such that 

tn = Ko{Sn)=n" 
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Proof We have go e~'^ implies 1 — Go{s) a{s)go{s) e"*^*) where ^ like V varies regularly with 
exponent 6. Write s„ = e*^" and t^'(s„) ^ e^'"^'^"^ where r is a function with r'{t) — > 1 and r"{t) — > 
as t — > oo, and r := 1/A. It has been shown in [2] (Equation (1.13)) that 

Ko{s) = t ~ ce"^'^) s ^ oo, Lp{s) = rg(*(s)) ~ t*(s), 

for some positive constant c. This gives logt„ = logii'o(sn) ^ e^*"^"^"^. Since log log t„ = n^~*logn + 
log log n has increments which go to zero, so does 6r{an), and hence also o-„ since r' tends to one. It 
follows that s„+i ^ s„. If 

Choose SnTO„-i = Sn-i- Then the cube [s„_ie, s„+ie] is a union of 2"^ blocks in the partition, and so 
is the cube [tn-ie, tn+ie] in z-space; e = (1, . . . , 1) denotes a vector of ones in M**. The union U of these 
latter cubes has the property that the scaled sets U/tn converge to (0,00)"^ for f„ — > 00 if t„ << tn+i- 

Theorem 3.7. Assume the standard set-up with the excess m,easure p 0/ the original distribution which 
does not charge the coordinate planes. Let p he an excess m,easure on \ {0} with marginal densities 
\/\t\^'^^ , and assume that for each orthant Qs, S G { — 1, 1}'', the restrictions of p and p to Qs have the 
same univariate marginals. One may choose F such that its marginals are tail asymptotic to Fq, such 
that the sample clouds converge to the Poisson point process N with mean measure p, and such that the 
sample clouds from the df G = F o K converge onto the limit set Ex^ in (2.8). 

Proof We sketch the construction. Choose F e L)\{p) with density / such that the sample clouds from 
F scaled by c„ converge to N. For 5 G {—1, l}*^ let Ug be the image of the union U in Qg by reflecting 
coordinates for which Sj = — 1 (see Figure 1 for an illustration). Let / agree with / on the 2"^ sets Ug and 
with / elsewhere, so that, by the remark above on the convergence of the scaled sets U, f and / differ 
only on an asymptotically negligible set. Alter / on a bounded set to make it a probability density. Then 
the sample clouds from F scaled by c„ converge to N. In the corresponding partition (A„) on x-space 
we only change the measure on the "tiny" blocks [sn-i, Sn+i]''' (with s„_i ^ Sn+i) around the positive 
diagonal, and their reflections. Hence the scaled sample clouds from G = F o K converge onto E\fi. ^ 

We now discuss the second case: the block partition (i3„) in z-space is regular, but {An) is not. As 
before, a block partition on x-space is determined by an increasing sequence of cubes s„C = [— s„, SnY, 
and for each n a symmetric partition of [— Sn,Sn] given by a finite sequence of points s„i < ••• < 
Snmn = in (0) Sn+i)- The image blocks are determined by f„ = i^o(sn) and tnj = Ko{snj)- This 
makes it convenient to define these quantities in terms of upper quantiles. Choose m„ = n, probabilities 
Pn = e~^, and write 

1 - Fo{t„) = 1 - Go(s„) =Pn 1 - Fo{tnj) = 1 - Go(Snj) = "Pn/j- (3.6) 

One may see tnj as a function T of j/npn, and so too Snj as a function S of \og{j/npn)- 

tnj = T{je^ln) Snj = S{y/n + \ogj - logn) j = 1, . . . , n, n = 1, 2, . . . 
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(a) z-space (b) x-space (regular) 

Figure 1: Panel (a): Schematic representation of the four sets Us = U„[(5it„_ie, ^2in+ie] in M? for 
6 = {51,62) € {—1, 1}^; see Theorem 3.7. Panel (b): the image of the sets Ug in x-space. The division 
points of the partitions (j4„) and (-B„) are s„ = n^logn and tn = e*" = n"^ respectively. Use 
Proposition 3.6 with Ko{s) = e^, s > so and e = 0.5. 

The increasing functions T and S vary regularly with exponents 1 /A and 1/9 since the inverse functions 
to 1 — -Fo and — log(l — Go) vary regularly in zero with exponents — 1/A and 1/6. It follows that 
T{{jn/n)e^)/T{e^) — > v}/^ if jn/n — > u e [0, 1], and hence tni/tn — » and the maximal increment 
tnj — tnj-i, j = 2, ...,n, is o{tn)- So the block partition in z-space is regular. However 
since logn = o{^/n) implies S{\/n — log n)/S'(\/n) — > 1- Figure 2 depicts sequences of cubes s„C 
and tnC in on which partitions and (-B„) are based in the special case when s„ = ^/n and 
tn = K(,{sn) = e*" = e^, along with subintervals [—e^/n,e^/n] x {e"} in z-space mapping onto 
[—y/n + logn, ^/n — logn] x {y/n} in x-space, which correspond to the partition blocks intersecting the 
coordinate axes. 

Theorem 3.8. Assume the standard set-up. There exists a df F such that the original df F and F have 
the same asymptotics, and the scaled sample clouds from the corresponding meta df G converge onto the 
diagonal cross Eqo in (3.3). 

Proof Let (-Bn) be the regular block partition in z-space as above. Construct a density / by deleting 
the mass of the original df F in the blocks i?„ which intersect one of the d coordinate planes, except 
for the block containing the origin, where we increase the density by a factor c > 1 to compensate for 
the loss of mass elsewhere. The new density / agrees with / on every block which does not intersect a 
coordinate plane. The relation i„i/t„ implies that / and / agree outside a vanishing neighborhood 
around the coordinate planes. On the other hand the relation s„i ^ s„ implies that g and g only agree 
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(a) z-space (regular) (b) x-space 

Figure 2: Panel (a): blocks [-t„,t„]'^ = [-e^,e^]^ n e {2, 4, 6, . . . , 100} in z-space with the sets 
[—e^/n,e^/n] x [e^"~-^, e"^] and their reflections, which correspond to the blocks of partition (B„) 
intersecting the axes, deleted. Panel (b): image of the blocks {[—^/n, and deleted blocks {[—\/n + 

logn, ^/n — logn] x [\/n — 1, \/n\) from Panel (b) in x-space under the meta transformation Ko{s) = e*, 
s > So- 

on a vanishing neighborhood around the diagonals. Thus f{rnWn) = f{rn^n) eventually for r„ — > oo if 
w„ converges to a vector w with non-zero components. Conversely, if r„ — > oo and w„ converges to a 
vector w which does not lie on one of the 2'^ diagonal rays, then 5(r„w„) = eventually. 

The function / satisfies the same limit relation f{rnWn)/f{rne) h{w) for r„ — > oo as /, provided 
w = lim„_>oo w„ has non-zero coordinates. Dominated convergence, by the inequality / < / outside tiC, 
gives L^-convergence outside centered balls. It follows that F and F have the same asymptotics. The 
scaled sample clouds from g converge onto the diagonal cross £^oo- If 

The incompatibility of the partitions (An) and = {K{An)) introduced in this section gives one 

technical explanation for the peculiar sensitivity of the limit shape for the meta distribution. If we regard 
the atoms of the partition (Bn) as nerve cells, then regularity of (j4„) will make the region aroimd the 
coordinate planes in z-space far more sensitive than the remainder of the space, and it is not surprising 
that cutting away these regions has drastic effects on the limit. 

3.2 Mixtures 

For a wide class of star-shaped sets A c [—1, 1]"^ it is possible to alter the original density so that the 
limit set E of the new meta density is the union A U Eqo, where Eqq is the diagonal cross in (3.3). 
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Theorem 3.9. Assume the standard set-up. Let A be a star-shaped closed subset of the unit cube [—1, 1]*^ 
with a continuous boundary and containing the origin as interior point, and let Eqq be the diagonal cross 
in (3.3). There exists a df F with the same asymptotics as F, such that the scaled sample clouds from 
the meta distribution G converge onto the set E = Au Eqo . 

Proof Let G = F o K where F has marginal densities /o, and let G° = F° o K where F° has continuous 
marginals F° with lighter tails than Fq, see (3.2). Let dF agree with d{F + F°) outside a bounded set. 
The sample clouds from F° scaled by c„ converge onto {0} since n(l — JJ(ec„)) + nJJ(— ec„) ^ as 
n — > oo, j = 1, . . . , d, e > 0. So F and F have the same asymptotics. 

Let UA denote the gauge function of A (or its interior), and let 5f°(x) = g'*(nA(x)) for a continuous 
decreasing positive function g^ on [0, oo). We assume that gf* varies rapidly. The function g° is continuous 
and < 5f°(x) < ^(x). Let dG agree with dG{x) + g°{x)dx outside a bounded set. We may assume that 
dG does not charge coordinate planes {xj = c} and charges all coordinate slices {ci < xj < C2}. Then 
the marginals are continuous and strictly increasing. 

Choosing F to have the same asymptotics as the original df F, and so that the scaled sample clouds 
from G converge onto the diagonal cross Eqq (see Theorem 3.8 and Proposition 3.1), we obtain F and F 
with the same asymptotics, and convergence of the scaled sample clouds from G onto Eqo U A. ^ 

Example 2. Let g be a density with cubic level sots: .g(x) = g*(||x||oo) with as in the above theorem. 
The marginal densities go are symmetric and equal, and asymptotic to (2s)''^^g, (,s). Let .go(fon) ^ l/n. 
We may choose g, such that fe„ ^ fe„ and ngo(feri) 0, (see Lemma A.l and Proposition A. 2 for details). 
It follows that the sample clouds from G scaled by 6„ converge onto the cube [—1,1]'', hence also the 
sample clouds scaled by 6„. 

In the above theorem, take F = F and G° = G. Then F has the same asymptotics as F and the 
scaled sample clouds from G = F o K converge onto the cube [—1, 1]''. 

If we choose the measure dF to be the image of the marginal dFo under the map 1 1— > te, then dF lives 
on the diagonal. The scaled sample clouds from F converge to the Poisson point process on the diagonal 
with intensity A/|t|'^+^ in the parametrization above, and the scaled sample clouds from G converge onto 
AuEoQ. If we choose A = Ex^g then the dfs G and G have the same asymptotics, but we may also choose 
A = E\*fi* for other values A* and 6* in (0, 00). 

We have shown that slight changes in F, changes which do not affect the asymptotics or the marginals, 
may yield a meta distribution G with the marginals of G but with different asymptotics. This makes it 
possible to start out with a Poisson point process N* and a star-shaped set E* , and construct dfs F* 
and G* with marginal densities /o and go such that the sample clouds from F* converge to A''* and those 
from G* converge onto E* . The only condition is that the mean measure p* of A^* is an excess measure 
with marginal densities A/ltl^"*"^ and that E* is a closed star-shaped subset of [—1, 1]'* containing the 2'^ 
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vertices and having a continuous boundary. 

Although the shape of the limit set is rather unstable under even slight perturbations of the original 
distribution, one may note the persistence of the diagonal cross as a subset of the limit set. Due to 
equality and symmetry of the marginals, clearly the points on the 2'^~^ diagonals in z-space are mapped 
into the points on the diagonals in x-space. However, Figure 1 shows that much larger subsets of the 
open orthants in z-space are mapped on the neighborhood of the diagonals in x-space. 

Proposition 3.10. Consider the standard set-up with the original df F having equal marginals. Let F 
he a df on M'' whose marginals are equal and symmetric and tail asymptotic to those of F. Assume that 
the sample clouds from F converge to a point process with mean measure p, where p charges (0,oo)'^. // 
sample clouds from the meta distribution G = F o K can he scaled to converge onto a limit set E then E 
contains the vertex e = (1, . . . , 1). 

Remark 2. Limit sets are always star-shaped (see Proposition 4.1 in [7]). Hence, if the limit set E 
contains e, it also contains the line segment E^q = Eqq n (0, oo)'' joining e and the origin. 

Proof To see this wc again make use of the block partitions of Section 3.1; in particular consider 
the situation sketched in Figure 1. Due to symmetry, it suffices to restrict attention to the positive 
orthant. Consider cubes := [s„ — Ma(s„),s„ + Ma(s„)]'^ centered at diagonal points ,s„e for some 
M > 0, where a{s) is the scale function of the marginal df Gq. Recall that a'(s) and hence 
a{s)/s ^ as s ^ 00. These cubes are asymptotically negligible as C^/sn ^ {e} for ,s„ — > oo. The 
corresponding cubes in z-space are centered at the diagonal points t„e with t„ — Ko{sn) and given by 
K{C^) = [Ko{s„ - Ma{snj),Ko{sn + Ma{sn))Y =: [tn-iJ^n+iV- The von Mises condition on 
1 — Go with scale function a(s) and regular variation of (1 — Fo)~^ in zero with exponent — 1/A give 

^^t_^^ Ms„ -Mais.)) ^ il-F,)-Hl G.js,. Majs.))) 

lim ( ^-Gojsn) yA L(l-Go(s„-Ma(s„))) 
«„™oo \l - Go(s„ - Ma{sn))J L{1 - Go(s„)) 

^ / e-^(l-Go(g„)) yA L(l-Go(g„)) _^/;, 

sn^^K 1-GoM ^ L(e-^(l-Go(s„))) 

for a slowly varying function L. Similarly, t„+i/t„ e^/^, and thus C^ /t„ [e"^/^, e^/'^]'*. Note 
that for M large, this limit constitutes a large subset of (0, oo)'' in z-space in that it will eventually, for 
M large enough, contain any compact subset of (0, oo)*^. 

The points in C^ are mapped by the meta transformation K onto the points in G^, thus preserving the 
mass P{Z e G^} = P{X e G^}. This shows that most of the mass on (0, 00)° in z-space is concentrated 
on just a neighborhood of Eqq in x-space, and the assumption on p ensures that the limit of the scaled 
sample clouds contain the vertex e. ^ 
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3.3 Extremes and high risk scenarios 

There are different ways in which one can look at muhivariate extremes. The focus of the paper so far 
has been on a global view of sample clouds. We have seen that the asymptotic behaviour of sample clouds 
is described by a Poisson point process for heavy-tailed distributions, and by the limit set for light-tailed 
distributions. We would now like to complement these global pictures by looking more closely at the 
edge of sample clouds. In particular, we discuss asymptotic behaviour of coordinatewise maxima and of 
exceedances over hyperplanes, termed high risk scenarios in [4]. 

3.3.1 Biregular pairtition and additional notation 

We first introduce a biregular block partition, a partition which is regular in both z-space and in x-space. 
For simplicity, we shall restrict attention to a bivariate situation. Recall the block partition defined in 
terms of quantiles in (3.6), which was regular in z-space, but not in x-space. We now refine it to make 
it biregular. A typical block intersecting the positive horizontal axis in x-space has the form 

[s„_i,S„] X [-S„l,S„l] Sn-1 Sn S„i. (3.7) 

It is an elongated thin vertical rectangle almost stretching from one diagonal arm to the other. We 
subdivide it in the vertical direction into 2n congruent rectangles by adding the division points 

-(n - l)s„i/n, . . . , -Sni/n, 0, s„i/n, . . . , (n - l)s„i/n. 

The corresponding partition (A„) in x-space is regular, and so is the partition (B„) with Bn = K{An) 
in z-space, since it is a refinement of a regular partition. 

We shall be interested in three sets C, D, and O which are disjoint and fiU They are defined in 
terms of the biregular partition we have introduced above. Roughly speaking, C consists of the blocks 
around the coordinate axes, D consists of the blocks around the diagonals, and O is the remainder. Due 
to the symmetry of the partition, we distinguish four subsets of C associated with the four halfaxes: 
Cn with {0} X (0,00), Ce with (0, 00) x {0}, and their negative counterparts Cs and Cw- Similarly, 
Dne is the restriction of D to (0,oo)^ and analogously for the other quadrants Dse, Dsw and Dnw 
clockwise. Here is the assignment of the blocks An (or _B„) to the sets C, D and O: In the set Ce we 
put the rectangles in (3.7), for n > no, or rather the 2n blocks into which these rectangles have been 
subdivided. Similarly, for Cat, C'w and Cs- In the diagonal sets D we put all blocks in the ring R„ 
between successive squares, which are determined by subdivision points ±.s„j with j > j„ ~ We 
claim that Dne asymptotically fills up the open positive quadrant in z-space. For this we have to show 
that tnj^/tn 0. By definition t„j is the upper g-quantile for the df Fq for q ~ ne^^/j. Let t„ be the 
upper quantile for p„ = and r„ the upper quantile for q„ — ^Jne^^ . We claim that r^jtu ^ 0. 

This follows by regular variation of the tail of Fq because Pn/ln 0. So in z-space the diagonal set D 
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(a) z-space (b) x-space 

Figure 3: Panel (a): Sequences of rectangles [i„_i,t„] x [t^, t„] filling region Ce and of cubes \Ui,tn]'^ 
filling region D^e, and their reflections in other quadrants with t„ = and i„ = {\y^/n)tn for 
n G {1,2,..., 20}. Panel (b): The images of the sets in (a) in x-space under the coordinatewise map 
t = Kois) = e-' for n £ {10, 15, ... , 200}. 

asymptotically fills up the whole plane apart from the coordinate axes; in x-space the coordinate set C 
fills up the whole plane apart from the two diagonal lines; and there is still a lot of space left for the set 
O (note that tni/tnj„ — + 0). See Figure 3 for illustration. 

The asymptotic behaviour of the probability distribution F is known once we know the probability 
masspn of the atomic blocks An in x-space (or of Bn = K{An) in z-space) (see Propositions 3.3 and 3.4). 
The specification of the probability masses p„ of the atoms is not a very efficient way of describing the 
original distribution and the meta distribution. Instead we shall describe the asymptotic behaviour of dF 
on the diagonal sector Dne, and the asymptotic behaviour of the meta distribution dG on the coordinate 
sector Ce- We first consider the two pure cases, where all mass lives on one of the sets C or D. 

3.3.2 The region D 

The situation on the region D is similar to that discussed in Theorem 3.8. In z-space, let F* be a df 
in the domain of an excess measure p on K.^ \ {0}, where p has marginal intensities X/t^^^, A > 0, and 
no mass on the axes. The measure p might be concentrated on one of the two diagonals. Since we are 
concentrating on the sector Dne we shall assume that p charges the positive quadrant. Delete the mass 
outside D, and compensate by adding some mass in a compact set. The marginals of the new df F are 
still tail asymptotic to _Fb- The measures ndF, scaled properly, converge to p weakly on the complement 
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of centered balls. The sample clouds and the coordinatewise maxima converge under the same scaling. 
The exponent measure of the max-stable limit distribution is the image p+ of the measure p under the 
map z = (a;, y) (a; V 0, y V 0) = z+. 

The following proposition describes the relation between the exponent measures for an original df and 
the associated meta df. It is in line with e.g. Proposition 5.10 in [10]. One can determine the max-stable 
limit for a meta df from the max-stable limit for the original df and the marginals of the meta distribution. 
We shall refer to this fact as the invariance principle (for coordinatewise extremes) since the max-stable 
limit distributions have the same copula. For a df in the domain of attraction of an extreme value 
limit law with exponent measure p, we use the notation F G DA{p). 

Proposition 3.11. Let F G DA{p'^) be a continuous df on whose exponent m,easure has m,arginal 
intensities all equal X/t^^^ , A > 0. Let Go be a continuous symmetric df on M. whose tails are asymptotic 
to a von Mises function. If G is the m,eta df based on F luith marginals equal to Go. then G G DA(cr^), 
where 0"+ and p~^ are related by a coordinatewise exponential transformation, p~^ = K{a~^) with equal 
components 

K:ut^w, Wi = Ko{ui)=e^'l^ i = 1,2. 

Proof The map K is the limit of coordinatewise transformations mapping normalized coordinatewise 
maxima from G into normalized coordinatewise maxima from F, hence it is also a coordinatewise trans- 
formation. The von Mises condition on Gq implies that Go is in the domain of attraction of the Gumbel 
limit law exp{— e~^}, and thus <t+ has standard exponential marginals. The relations 

pt[t, oo) = cr+[s, oo) <^ f-^ = e"** <^ t = Kois) = e"'^, i = 1,2, 

for the marginals determine K. ^ 

In x-spacc Dj^e is a thin strip along the diagonal. However, it follows from Proposition 3.11, that 
the coordinatewise maxima in x-space, centered and scaled, converge to a max-stable distribution whose 
exponent measure cr"*" is the coordinatewise logarithmic transform of The measure (T+ lives on 

[—00,00)^ \ {(—00,-00)}. The restriction a to is the image of the restriction of p+ (or of p) to 

(0,«3)2. 

Let us now look at the asymptotic behaviour of the high risk scenarios and sample clouds for F and 
G. Let Hn = > c„} be halfplanes with direction ^„ = (a„,6„) of norm one, and c„ 00. For the 
heavy-tailed df e f (p), the domain of attraction of the excess measure p, the situation is simple. 

Proposition 3.12. If ^ ^ o,nd c„ — > 00 then Z^" W, where W has distribution dp^ = Indp/pH 
forH = {^>l}. 

Proof The map (u, v) 1-^ {x, y) — (c„m, c„7;) maps J„ = > 1} onto By assumption the probability 
measures dF scaled by c„ and multiplied by a suitable factor converges to p weakly on the complement 
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of centered balls, hence on H. This gives the result when = ^ for all n. For the general case, — * ^ 
use the continuity theorem [4], Proposition 5.13. ^ 

The same result holds for the light-tailed meta distribution G provided the limit direction ^ does not 
lie on one of the axes. 

Proposition 3.13. //Cn ^ C ^ (0, oo)^ and c„ ^ oo then a^^CX.^"^) = (X"^" — (6„, bn))/an => U where 
(bn, bn) € dHn and a„ = a(fe„) for the scale function a associated with the marginal density go. The limit 
U has distribution du^ = luda/uH for = > 0}. 

Proof This follows from the weak convergence of d-K = dG normalized by and multiplied by l/p„: 

'^n^M/Pn cr^/ao weakly on \ [— oo, c]^ c e M 

for pn = 7r(IR2 \ [6„,oo)2) and ao = a+{[-oo,oof \ [-oo,0]^). 1 

Convergence of high risk scenarios implies convergence of sample clouds with the same normalizations 
for halfplanes that satisfy P{Z e Hn} ~ 1/n as in [4], Section 14. For the heavy-tailed df F the 
convergence of the sample clouds Nn N weakly on > c} for all c > is no surprise since Nn N 
holds weakly on the complement of centered disks. For the light-tailed df G weak convergence on > c} 
for ^ e (0,00)^ follows from weak convergence on [— oo,cxd)^ \ [— c»,c]^. Convergence for horizontal 
halfspaces i?„ = {y > c„} presents a different picture, since we do not allow mass to drift off to the 
vertical line in —00. 

Proposition 3.14. Set Ht = {y > t} and at{u,v) = {tu,t + a{t)v). Then aj~^(X^«) {U,V) where 
U and V are independent, V is standard exponential, and U assumes only two values, F{U = —1} = 
p- = c-/c and F{U = 1} = p+ = c+/c. Here c- = <t+({-oo} x (0,oo)) = p+({0} x (l,oo)) = 
p((-oo, 0) X (1,00)), c+ = a-(M X (0,00)) = p((0,oo) x (l,oo)), and c = p(R x (l,oo)) since by 
assumption p does not charge the axes. If P{X e -ffr„} ~ ^/n then the sample clouds converge 
Mn = . . . , q:^^(X„)} M weakly on {v > c}, c G R. The Poisson point process M lives on 

two vertical lines, on {u = —1} with intensity p-e~'" and on {u = 1} with intensity p+e"". 

Proof It suffices to prove the second relation. The limit point process M with mean measure a above the 
line {v = —C} for C > 1 corresponds to the sample points Xi, . . . , X„ above the line {y = rn ~ Ca(r„)}. 
The corresponding points scaled by r„ converge to the vertex (1,1) of the limit set E. Under the 
normalization a^^ the horizontal coordinate converges to 1 , and hence the whole sample cloud converges 
to the projection of M onto the line {m = 1}. A similar argument holds for the points of the point 
process with mean measure a on associated with restriction of p to the quadrant (— oo,0) x (0, 00). 
These yield the points on the vertical line through (—1,0). Since p does not charge the vertical axis, this 
accounts for all points. ^ 

On the region D there is a simple relation between the Poisson point processes associated with the 
exponent measures and the Poisson point processes associated with the high risk scenarios. 
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3.3.3 The region C 

Wc shall consider the counterpart of the class J-'\ for light-tailed densities g on the plane. The density g is 
assumed to be imimodal, continuous and to have level sets {g > c} which are scaled copies of a bounded 
open star-shaped set S in the plane with continuous boundary. Such densities are called homothetic. 
They have the form g — g*{ns) where is a continuous decreasing function on [0,oo), the density 
generator, and ns is the gauge function of the set S (sec Table A. 3). Wc shall assume that the density 
generator is asymptotic to a von Mises function, and that 5 is a subset of the square (—1, 1)^. 

If S is the open unit disk, or more generally a rotund set (i.e., convex with a boundary having 
positive curvature in each point), then the high risk limit scenarios exist for all directions, and are Gauss- 
exponential. The associated sample clouds in these directions may be normalized to converge to a Poisson 
point process with Gauss-exponential intensity e~"^/^e~^/\/27r on ]S? for appropriate coordinates u,v. 
See Sections 9-11 in [4]. 

If the set 5 is a convex polygon with one boundary point q = {qi, 1) on the line {y = 1} with \qi\ < 1, 
then the horizontal high risk scenarios converge. The associated excess measure has density h, where h 
has conic level sets {h > e~*} = C + tci for an open cone C in the lower halfplane. The cone C describes 
the asymptotic behaviour of the set 5 — q in the origin. 

If 5 is the square (—1, 1) then the horizontal high risk scenarios, normalized by at{u, v) = {tu,t+a{t)v) 
converge. The associated excess measure has density l[_i_i](u)e~'', and vanishes outside the vertical strip 
[-1,1] X M. 

Our density g is defined in terms of a bounded star-shaped set S and the density generator g» which 
describes the behaviour of g along rays up to a scale factor depending on the direction. On the other hand 
we want light-tailed densities with prescribed marginals equal to go- In general it is not clear whether 
for a given set S there exists a density generator (7* which produces the marginals go- 

Example 3. Let go be the standard normal density. Let g = g^{ns)- Given a bounded open convex 
set S can one find a density generator such that g has marginals 50? If 'S is the unit disk, or an 
ellipse symmetric around the diagonal, then g*{r) = ce~'^^^^ will do. If S is the square (—1,1)^, then 
one may choose g*{r) ce~'"^/^/r so that the tails of the marginals agree with go, see Section A. 2 in [2]. 
By altering 5 on a square [— M, M]^ one may achieve standard normal marginals. If S is the diamond 
spanned by the unit vectors on the four halfaxes, then one can choose 5*(r) cre~^^/^ so that the tails 
of the marginals agree with go, see [9]. These three density generators are different. It is not clear how 
to combine the asymptotic behaviour of these examples. 

Let us take a mixture of these three densities with weights 1/3 each. Set at{u, v) = {tu, t+a{t)v) where 
a{t) = 1/t is the scale function associated with the normal marginal go- Since under the normalization 
at all three distributions have high risk limit scenarios {U, V) with U and V independent and V standard 
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exponential, and since the normalizations for the vertical coordinate V are determined by the vertical 
marginal go which is the same in the three cases, we concentrate on the horizontal component. For the 
density with square level sets, U is uniformly distributed on [—1, 1], whereas in the other two cases U has 
a point mass at the origin. Hence, one obtains for the mixture a uniform distribution on [—1, 1] with a 
point mass of weight 2/3 in the origin. If we use the partial compactification of the plane in [6] then it is 
also possible to obtain limit distributions on [— oo, oo] for the horizontal component. There are several: 
A centered Gaussian density with an atom of weight 1/3 in the origin and two atoms of weight 1/6 in 
±oo; a Laplace density with two atoms of weight 1/3 in ±oo; an atom of weight 1/3 or 2/3 in the origin 
with the remaining mass fairly divided over the two points in oo; all mass in the points in oo. 

Now consider an open ellipse E, symmetric around the diagonal, which has a unique boundary point 
on the line {x2 = 1} in (p, 1) for some p £ (0, 1). Its reflection E' in the vertical axis has the boundary 
point {—p, 1). The union S = EU E' is a star-shaped set. What does the excess measure for horizontal 
halfspaces look like? If we zoom in on {p, 1) we obtain a Gauss-exponential measure, but the measure 
around the point {—p, 1) moves off to — oo. If we want weak convergence on horizontal halfspaces we have 
to use the normalization at above. The limit measure now lives on the two vertical lines u = ±p and has 
the same exponential density e"" on each by symmetry. 

If 5 C (—1,1)^, and the boundary of the star-shaped set S contains points on the interior of the 
four sides of the square, but does not contain any of the vertices, then the components Xi and X2 of 
the vector X with density g are asymptotically independent, and so are the pairs (—^1, X2), (Xi, —X2), 
(— Xi,— X2), see [3]. This also holds if S is the whole square. So if g is the meta density based on a 
heavy-tailed density in I'(p), then the excess measure p lives on the four halfaxes, and the distribution 
may be restricted to the region C. 

If the high risk scenarios for horizontal halfplanes converge with the normalizations at{u, v) = {tu, t + 
a{t)v) of Proposition 3.14 then the limit vector {U, V) has independent components and U lives on the 
linear set Ei determined by the intersection of the boundary of the limit set E with the horizontal line 
{v = 1}. Note that in this setting the limit set E is the closure of S, see [3]. 

Proposition 3.15. Let G be a bivariate df with marginal tails all asymptotic to the von Mises function 
with scale function a. Let X* denote the high risk scenario X^ for the horizontal halfplane H = 
{y > t}. Let at{u, v) = {x, y) = {tu, t + a{t)v). Suppose q:7^(X*) U = V) for t — > 00. Then U and 
V are independent, V is standard exponential and \U\ < 1. Let tjj{bn) = logn. If the sample clouds from 
G scaled by bn converge onto the compact set E then F{{U, 1) G E} = 1. 

Proof The distribution of V is determined by the marginal df G2 and is exponential since 1 — G2 is 
asymptotic to a von Mises function, see [4], Proposition 14.1. Light tails of the marginal Gi imply 
(1 - Gi{eH))/{l - Gi{t)) for any e > 0, and hence F{U > e"} = 0. Similarly for the left tail. 
Let a denote the excess measure extending the distribution of {U,V) to M^; see [4], Section 14.6. Since 
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ct^^ at_^_sa,{t) — * Ps where ps{u,v) = {u,v + s), it follows that a{A + {0,t)) = e~*cr{A) for all Borel 
sets A in M?. See Proposition 14.4 in [4]. This implies that cr is a product measure, and hence that 
U and V are independent. Let (3n{u,v) = {bnU,bn + bnv). The sample clouds {/3~^(Xi), . . . ,/3~^(X„)} 
converge onto the translated set E— (0, 1). Now apply the additional normalization 7„('U, v) = {u, anv/bn) 
where a„ = a(6„) implies Un/bn 0. Then ')n^Pn^ = (x~^ where an{u,v) = {bnU,bn + a„w). Hence 
the renormalized sample clouds converge to the Poisson point process with mean measures a, and the 
restriction of a to {v > 0} lives on Eq x [0, oo) where Eq is the set {u \ {u, 1) e E}. ^ 

3.3.4 The combined situation 

In a mimbor of examples above of distributions on C the high risk scenarios for horizontal halfplanes, 
normalized by at : {u, v) (.x, y) = (tu, t + a{t)v), converge to a random vector {U, V) with independent 
components, V is standard exponential and U lives on the interval [—1, 1]. Recall from Proposition 3.14 
that for meta distributions on D these high risk scenarios with the same normalization also converge to 
a random vector with independent components. The vertical component is again standard exponential; 
the horizontal component lives on the two point set {—1,1}. The associated excess measmes are product 
measures on the vertical strip [—1, 1] x M. The density along the vertical line is cqc^^' /c. Now consider 
the sum of these two light-tailed distributions. We can alter this sum on a compact set to make it into a 
probability measure. The excess measure a associated with the high risk limit distribution for this new 
light-tailed distribution has the same structure: It is a product measure on the vertical strip [—1,1] x M, 
the vertical component has density ce~^, the horizontal component is a probability measure a* on [—1, 1] 
with atoms p_ and in the points ±1. 

If we look at the heavy-tailed distributions we see that the excess measures p for distributions on D 
live on the complement of the axes, and for distributions on C live on the axes. So here the sum has 
an excess measure which has mass both on the axes, and on the complement. The restriction of this 
measure to the set above the horizontal axis is characterized by the projection on the vertical coordinate 
with density cX/r^'^^ and a probability measure p* on the horizontal line, the spectral measure (see [4], 
Section 14.8) with the property that p*[c,oo) = p{E)/c, where E = {{u,v) \ u > cv.v > 1} is the set 
above the horizontal line {v = 1} and to the right of the ray through the point (c, 1). The probability 
measure p* has an atom of weight po in the origin. 

Proposition 3.16. The weights p-,p+,po above have sump- +po > 1- 

Proof If the horizontal component U on [—1,1] from the high risk limit scenario due to the light- 
tailed density on C lives on the open interval (—1,1) then p„ -f po + P+ — 1; ^nd pi = Ci/c as in 
Proposition 3.14 where now c_ = /?((— oo, 0) x [1, oo)), cq = p({0} x [1, oo)), c+ = p((0, oo) x [1, oo)), and 
c = c_ -l-co -I- c+ = p(M X [l,oo)). t 

There is a simple relation between high risk scenarios for horizontal halfspaces for the heavy-tailed 
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density and the light-tailed meta density on in terms of the meta transformation provided the halfs- 
paces have the same probability mass. The inverse maps such a z-halfspace H into an x-halfspace 
H'. In the limit for the excess measures this yields a coordinatewise mapping from z-space to x-space: 

(21,2:2) 1-^ {xi,X2) = (sign(2i),log22)- (3.8) 

The transformation for the horizontal coordinate is degenerate. Define the curve T as the graph of a; 1— > 
sign(a;) augmented with the vertical segment at the discontinuity in zero. There is a unique probability 
measure on F whose horizontal projection is p* and whose vertical projection is a*. This follows from 
the inequality in Proposition 3.16. The excess measures p on (0, 00)^ and a on are linked by the 
coordinatewise exponential transformation in Proposition 3.11. The excess measures p on M x (0,oo) 
and a on [—1,1] x M are linked by the coordinatewise map (3.8). The high risk limit scenarios are 
limits of vectors and where H' and H are corresponding horizontal halfspaces (with the same 
mass). These vectors are linked by the coordinatewise monotone transformation K restricted to H'. In 
such a situation only a very limited class of transformations is possible in the limit, apart from afiine 
transformations only power transformations, the exponential and its inverse the logarithmic, and four 
degenerate transformations amongst which the transformation sign and its inverse, see [1] Chapter 1 for 
details. 

4 Discussion 

In situations where chance plays a role the asymptotic description often consists of two parts, a deter- 
ministic term, catching the main eSect, and a stochastic term, describing the random fluctuations around 
the deterministic part. Thus the average of the flrst n observations converges to the expectation; under 
additional assumptions the diff'erence between the average and the expectation, blown up by a factor 
^/n, is asymptotically normal. Empirical dfs converge to the true df; the fluctuations are modeled by a 
time-changed Brownian bridge. For a positive random variable, the n-point sample clouds Nn scaled by 
the 1 — 1/n quantile converge onto the interval [0, 1] if the tail of the df is rapidly varying; if the tail is 
asymptotic to a von Mises function then there is a limiting Poisson point process with intensity e~*. 

Convergence to the first order deterministic term in these situations is a much more robust affair than 
convergence of the random fiuctuations around this term. So it is surprising that for meta distributions 
perturbations of the original distribution which do not affect the second order fiuctuations of the sample 
cloud at the vertices may drastically alter the shape of the limit set, the first order term. This paper tries 
to cast some light on the sensitivity of the meta distribution and the limit set E to small perturbations 
of the original distribution. 

Bivariate asymptotics are well expressed in terms of polar coordinates. Two points far off are close 
together if the angular parts are close and if the quotient of the radial parts is close to one. This geometry 
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is respected by certain partitions. A partition is regular if points in the same atom are uniformly close as 
one moves out to infinity. Call probability distributions equivalent if they give the same or asymptotically 
the same weight to the atoms of a regular partition. Equivalent distributions have the same asymptotic 
behaviour with respect to scaling. 

This paper compares the asymptotic behaviour of a heavy-tailed bivariate density with the asymptotic 
behaviour of the meta density with light-tailed marginals. Small changes in the heavy-tailed density, 
changes which have no influence on its asymptotic behaviour, may lead to signiflcant changes in the 
asymptotic behaviour of the meta distribution. We show that regular partitions for the heavy-tailed 
distribution and for the light-tailed meta distribution are incommensurate. The atoms at the diagonals in 
the light-tailed distribution fill up the quadrants for the heavy-tailed distribution; atoms at the axes in the 
heavy-tailed distribution fill up the four segments between the diagonals for the light-tailed distributions. 
Section 3.1 shows how equivalent distributions in the one world give rise to different asymptotic behaviour 
in the other. 

In our approach the asymptotic behaviour in both worlds is investigated by rcscaling. In the heavy- 
tailed world one obtains a limiting Poisson point process whose mean measiire is an excess measure p which 
is finite outside centered disks in the plane; in the light-tailed world the sample clouds converge onto a 
star-shaped limit set E. The only relation between p and E is the parameter A. This parameter describes 
the rate of decrease of the heavy-tailed marginal distributions; it also is one of the two parameters which 
determine the shape of the limit set E. The measure p describes the asymptotics for extreme order 
statistics; the set E for the intermediate ones. In Section 3.2 it is shown that it is possible to manipulate 
the shape of the limit set E without affecting the distribution of the extremes. 

There are two worlds, the heavy-tailed and the light-tailed; the bridge linking these worlds is the meta 
transformation that (in our approach) maps heavy-tailed distributions into light-tailed distributions. 
Coordinatewise rmiltivariate extreme value theory with its concepts of max-stable laws and exponent 
measures is able to cross the bridge, and to describe the asymptotic theory of the two worlds in a 
unified way. The exponent measures of the heavy-tailed world are linked to the light-tailed world by a 
coordinatewise exponential transformation. 

A closer look reveals a different universe. The two worlds exist side by side like the two sides of a sheet 
of paper, each having its own picture. On the one side we see a landscape; on the other the portrait of a 
youth. Closer inspection shows an avenue leading up to a mansion in the far distance in the landscape, 
and a youth looking out from one of the windows; in the black pupils of his eyes on the reverse side of the 
paper we see a reflection of the landscape. In Section 3.3 it is shown that it is possible to disentangle the 
heavy- and the light-tailed parts of a distribution by using the marginals to deflne biregular partitions. 

The heavy-tailed distribution gives information about the coordinatewise extremes. The light-tailed 
distribution gives this information and more: a limit shape, and limit distributions for horizontal and 
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vertical high risk scenarios. The Umit shape is a compact star-shaped subset of the unit square; the 
horizontal and vertical high risk limit scenarios have independent components {U,V). In appropriate 
coordinates, V is standard exponential and U is distributed over the interval [—1, 1]. 

Biregular partitions make it possible to combine a representative of the heavy-tailed distribution (with 
no excess mass on the axes) with a representative of the light-tailed distribution (with asymptotically 
independent components) into one probability distribution. The excess measure of the combined distri- 
bution is the sum of the two individual excess measures; the limit set is the union of the two individual 
limit sets. The biregular partition is finer than the regular partition for the light-tailed distribution. It 
not only sees the limit set, it is fine enough to discern the high risk asymptotics in all directions. 

While the asymptotics of heavy-tailed multivariate distributions is well understood and nicely reflected 
in the asymptotic behaviour of the copula at the vertices of the imit cube, for light-tailed distributions 
there are still many open areas. The possible high risk limit scenarios for light-tailed distributions have 
been described; see [4], Section 14. However, their domains of attraction are still imexplored territory. 
This paper shows that the asymptotics of heavy-tailed multivariate distributions do not suffice to describe 
the asymptotic behaviour of the light-tailed meta distributions. 

One possible explanation for the loss of information in crossing the bridge between heavy and light tails 
is the highly non-linear nature of the meta transformation. Coordinate hyperplanes are preserved, and so 
are centered coordinate cubes because of the equal and symmetric marginals. Geometric concepts such as 
sphere, convex set, halfspace and direction are lost. The meta transformation maps rays in the x-world 
which do not lie in one of the diagonal hyperplanes {x, = Xj } into curves which are asymptotic to the 
halfaxis in the center of the sector bounded by the diagonal planes. Similarly, the inverse transformation 
maps rays which do not lie in one of the coordinate planes into curves converging to the diagonal ray in 
the center of the orthant bounded by the coordinate planes. In the bivariate setting there is a clear duality 
between diagonals and axes. Heavy-tailed asymptotic dependence reduces to co- or counter- monotonicity 
in the light-tailed world; horizontal and vertical high risk limit scenarios in the light-tailed world reduce 
to asymptotic independence in the heavy-tailed world. Biregular partitions allow us to see the combined 
effect, by describing the two worlds side by side. 

The asymptotics of light-tailed densities with given marginals is not well understood. We hope to 
return to this topic later. Similarly it is not clear what role is played by the remainder set O in our 
decomposition CD DuO. 
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A Appendix 

A.l Supplementary results 

Lemma A.l. Let g be a positive continuous symmetric density which is asymptotic to a von Mises 
function . There exists a continuous unimodal symmetric density gi such that for all c G (0, 1) 

gi{s)/g{s)^0 gi{cs) / g{s) ^ oo s ^ oo. 

Proof Let M„(s) = tp{s) - ipis - s/n), and let M*{s) = mmt>s Mnit) for n > 2. Each function M* 
is increasing, continuous and unbounded (since t/a{t) oo), and for each s > the sequence M*(s) is 
decreasing. There exists a continuous increasing unbounded function h such that 

Hm M„(s) - 6(s) = oo n = l,2, .... (A.l) 

Indeed, define M*{s) = M*{s) on [a„,6„] = {M* e + 1]}, and M*{s) = n on [6„_i,a„]. Then 
M*{s) is increasing and M*(s) < M*{s) eventually for each n > 2. Set b{s) = M*(s)/2 to obtain (A.l). 
The function gi = e~^^ with ipi{s) = ^j{s) + b{s) is decreasing on [0,oo) and continuous, and 6(s) — > oo 
implies gi{s)/g{s) ^ for s ^ oo. For c = 1 — l/m the relations 

ip{s) — tpi{cs) = ip{s) — tp{cs) — b{cs) > Mm{s) — b{s) — > oo 

hold, and yield the desired result. ^ 

Proposition A. 2. Let g^ be a continuous positive symmetric density which is asymptotic to a von Mises 
function e~'^. Choose r„ such that gd{s)ds 1/n. Let gi be the probability density in Lemma A.l. 
There exists a unimodal density g{x) = gi*(||x||cxD) on R"^ with cubic level sets and marginals gi. The 
sample clouds from the density g, scaled by rn converge onto the standard cube [— l,!]"^. The functions 
/i„(u) = nrf^g{rnu) are unimodal with cubic level sets. They satisfy 



hn{u) 



oo ue(-i,i)'' 
u^[-l,l]''. 



Let E be a closed subset of C = [—1, 1]"^, containing the origin as interior point, star-shaped with contin- 
uous boundary. Set ce = |-B|/2'^. Then 5b(x) = g*{nE{^) / ce is a probability density, and the sample 
clouds from gE scaled by rn converge onto the set E. 

Proof Existence of g follows from Proposition A. 3 in [2]. Let Gi be the df with density gi. Then 
n(l — Gi(cr„)) for c > 1, and n(l — Gi(cr„)) oo for c G (0, 1). Let tt be the probability distribution 
with the unimodal density g. The limit relations on the marginal dfs Gi imply that mr{Bn) oo for 
the block i?„ = [— 2r„, 2r„]'^~^ x [cr„, 2r„] for any c € (0, 1). Since ft,„ is unimodal with cubic level sets it 
follows that /i„ — > oo uniformly on [—a, c]'^ for any c e (0, 1). (Since hn{cl) < k implies nir^Bn) < (2c)'^fc.) 
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The area of a horizontal slice of the density gs sd, level u/ce > is less than the area of the horizontal 
slice of g at level y, but the height of the slice is proportionally more by the factor ce- So the slices have 
the same volume. The level sets of the scaled densities are related: 

{h-E > t/cE} = rE {h>t} = rC. 

So the function He mimics the behaviour of h = he- If 

A. 2 Summciry of notation 

This section is intended to provide a guide to the notation. Throughout the paper, it is convenient to 
keep in mind two spaces: z-space on which the heavy-tailed dfs F,F*,... are defined, and x-space on 
which the light-tailed dfs G,G*,... are defined. Table A.l compares notation used for mathematical 
objects on these two spaces. Table A. 2 can be consulted while reading Section 3 in order to keep track of 
various symbols used to distinguish original and meta distributions with certain properties and purpose. 

Table A.l: Symbols used to distinguish various objects of interest in z-space and in x-space. Notation 
for marginals assumes that all marginal densities are equal and symmetric as in (3.1). 



z-space x-space Comments 



F and / 


G = FoK andg 


joint df and density 


Fq and /o 


Go= Fqo Kq and go 


marginal df and density 


TT 




probability measures 






mean measures of scaled sample clouds 


Z, Zi, Z2, . . . 


X, Xi,X2, . . . 


random vectors 


Nn 


Mn 


scaled n-point sample clouds 


N: Poisson point process 


E: limit set 


limit of scaled sample clouds 


c„ : 1 - i=o(c„) ^ 1/n 


bn ■■ - log(l - Go{bn)) ^ logn 


scaling constants 


{Bn), tn = Ko{Sn) 




block partitions, division points 
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Table A. 2: Symbols used in Section 3 to discuss robustness and sensitivity properties of meta distributions. 



Where not explicitly mentioned, notation for the associated meta df is analogous, e.g. G = 


similarly for other objects /, E, etc. 


Symbols 


Description 


F,G 


dfs satisfying the assumptions of the standard set-up 


F*, G* 


dfs constructed to have properties (P.1)-(P.4) and yield the same asymptotics 




for F and F* , or G and G*; 


F 


any original df with marginals F^ 


F 


df obtained by changing F to have marginals Fj tail equivalent to Fq 




(e.g. mixtures in Section 3.2) 


F" 


original df with lighter marginals than F; see (3.2) 



Table A. 3: Miscellaneous symbols. 

Syu)l)()l Doscriptiou 

f >^ f ratios /(x)//(x) and /(x)//(x) are bounded eventually for ||x|| — > oo 
/ ~ / K^)/f{^) ^ 1 for ||x|| ^ oo 

ud gauge function of set D: n£)(x) > for x 7^ 0, 2? = {ud < 1} and n£)(cx) = cn£)(x) 

for X e M**, c > 
e a vector of ones in 

D tJic 01)011 Euoli(l(vui uiiil 1)m]1 in 3:'' 
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